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PACS 7 5 . 4 5 . + j - Macroscopic quantum phenomena in magnetic systems 
PACS 7 5 . 5 . Xx - Molecular magnets 

PACS 7 6 . 3 . -V - Electron paramagnetic resonance and relaxation 

Abstract. - We measure magnetization changes in a single crystal of the single-molecule magnet Fes when 
exposed to intense, short (<20 ^s) pulses of microwave radiation resonant with the m = 10 to 9 transition. 
We find that radiation induces a phonon bottleneck in the system with a time scale of ~5 /is. The phonon 
bottleneck, in turn, drives the spin dynamics, allowing observation of thermally assisted resonant tunneling 
between spin states at the 100-ns time scale. Detailed numerical simulations quantitatively reproduce the data 
and yield a spin-phonon relaxation time of Ti ~40 ns. 



> 

m 

o 
o 



X 



Since the discovery of resonant tunneling between spin 
states more than a decade ago [1,2], single-molecule magnets 
(SMMs) have been intensively studied. In the past few years, 
much effort has focused on the behavior of SMMs in the pres- 
ence of microwave/millimeter-wave radiation as a way to un- 
derstand the fundamental spin dynamics and the coupling of 
the spin to its environment [3-12]. This work has been moti- 
vated, in part, by the possibility that these systems could serve 
, as qubits [13]. 

At low temperatures (<10 K), the Feg SMM behaves as a 
spin-10 object with uniaxial anisotropy. The spin dynamics can 
be well described by the Hamiltonian 



n 



-DSl + E{Sl - Sl) + C{SX + St)- gfieS ■ H, (1) 



where the anisotropy constants D, E, and C are 0.292 K, 0.046 
K, and -2.9 x lO^^K, respectively, and 5 = 2 [14-16]. The 
first term produces a double-well potential for the spin's orien- 
tation, making the "up" and "down" directions (relative to the z 
axis) lowest in energy and producing a ~ 25 K baiTier between 
the two orientations [17]. A magnetic field H along the z axis 
makes one well lower in energy than the other, as illustrated 
in the insets to Fig. 4. There are 25h-1 = 21 energy levels for 
this S* =10 system. The second and third terms in Eq. 1 break 
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its cylindrical rotational symmetry and result in tunneling be- 
tween levels. Resonant tunneling occurs when the magnetic 
field causes levels in opposite wells to align. 

Recent studies of the magnetization dynamics in a radia- 
tion field have attempted to observe radiation-induced dynam- 
ics (such as Rabi oscillations) that would allow a direct deter- 
mination of the lifetimes of excited spin states and the dephas- 
ing time T2. In previous work on Feg, we showed that such 
efforts are complicated by the fact that resonant radiation heats 
the sample and drives the spins and lattice out of equilibrium 
on the millisecond timescale [6], a process that can be quanti- 
tatively described [8]. To circumvent this heating phenomenon, 
we have done experiments at a much shorter time scale in which 
an Feg sample is subjected to intense, short microwave pulses. 
Such a time-domain technique allows us to investigate the mag- 
netization dynamics as a temporal sequence of transitions be- 
tween spin levels. We observe the development of a phonon 
bottleneck that limits the spin-phonon relaxation and plays an 
essential role in the magnetization dynamics. The phonon bot- 
tleneck drives subsequent relaxation and allows us to directly 
observe the thermally assisted resonant tunneling process at 
time scales of <100 ns. Through detailed numerical simula- 
tions, we deduce an excited-state lifetime (Ti) of ^40 ns for 
this system. 

A single crystal of Feg was mounted in a cylindrical reso- 
nant cavity with Q^6200. The radiation field Hi at the sam- 
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pie position was estimated to be 1 Oe. The sample was fixed 
above a lithographically defined inductive pick-up loop, with 
the crystal's b axis parallel to the plane of the loop to maximize 
flux coupling. The emf induced in the loop when the sample's 
magnetization changed was measured with a SQUID voltmeter 
coupled to a room-temperature amplifier. Further details of the 
experimental setup are described elsewhere [7]. To character- 
ize our sample, we performed standard electron-spin resonance 
(ESR) reflection spectroscopy (not shown), which allowed us 
to determine the sample's orientation angles 9 = 37.5° and 
= 138° (defined hy S ■ H = H{S, cose + 5*:^ sin 6* cos + 
Sy sin 9 sin ip)), consistent with the directly measured orienta- 
tion of the crystal. 

Figure 1 shows the induced emf produced by our sample 
when a 20-/is pulse of 117.566-GHz radiation is applied. At 
a magnetic field of 1800 Oe, where the radiation is resonant 
with the transition between states m = 10 and m = 9, we ob- 
serve a clear signal (proportional to dM/dt), punctuated by sud- 
den jumps at the times when the radiation is turned on and off. 
In contrast, data taken at zero field, where the radiation does 
not couple to the sample, shows no detectable signal. For ref- 
erence, the inset of Fig. 1 shows the 1800-Oe data after nu- 
merical integration to obtain AM as a function of time. The 
data indicate that during the radiation pulse the magnetization 
decreases, as population is pumped out of the m = 10 state. 
After the radiation is turned off, the magnetization continues 
a downward trend, an effect we have previously characterized 
as being due to the spins and lattice having been driven out 
of equilibrium by the radiation [6, 8]. On top of this overall 
downward trend, there is a decay immediately after the radia- 
tion is turned on or off. We fit the data after the radiation is 
turned off to an exponential decay plus a constant term (the 
latter to account for the slow heating that occurs with a time 
scale of ^1 ms [6, 8]). The results of the fit are shown by the 
dashed line in the figure and yield a time constant of 4.4(3) 
/xs. We typically find a decay time of r^5 /is for this sample, 
with no systematic dependence on temperature, field or radia- 
tion power We interpret this relaxation as the signature of a 
phonon bottleneck [18, 19] in which emission of phonons dur- 
ing the decay from m = 9 to 10 leads to a build up in the pop- 
ulation of phonons resonant with that transition. The timescale 
for spins to populate the excited state during the radiation pulse 
or return to the ground state after the pulse is determined by 
the time for the phonon distribution to build up or decay away, 
respectively, the phonon-bottleneck lifetime Tpb. Evidence for 
phonon-bottleneck effects have been seen in other SMMs [20] 
and the effect has been suggested to occur in Feg [10]. Ours is 
the first direct measurement of Tpf, in a SMM system. Tpf, rep- 
resents the time scale either for the resonant phonons to ther- 
malize via nonlinear processes within the crystal or to escape 
into the environment of the cryostat. 

In order to study spin dynamics before the complication of 
lattice heating sets in, we focus on shorter pulses of 2-/is dura- 
tion. Some examples of the dM/dt signal as a function of time 
for different values of magnetic field at 2. 1 K (3.3 K) are shown 
in Fig. 2a (2b). To elucidate the dynamics, we analyze the data 
by taking several time slices of the curves in Fig. 2 (as well as 
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Fig. 1: (Color online) Emf (cx dM/dt) induced in our pick-up loop as 
a function of time when a 20-/is radiation pulse of frequency 1 17.566 
GHz is applied to an Feg crystal at 2.1 K in zero field and in a field 
of 1800 Oe, as designated. The dashed line shows an exponential fit 
to the 1800-Oe data after the end of the radiation pulse. The inset 
shows the change in magnetization AM obtained from numerically 
integrating the 1800-Oe data. 

similar curves at other field values, not shown) to measure the 
dM/dt "lineshape" as a function of field. The results of such an 
analysis are shown in Fig. 3. Figure 3a shows the data at 2.1 
K. At this temperature, the dynamics are essentially restricted 
to the m = 10 and 9 levels (see Fig. 4a inset). When the radia- 
tion is on (solid symbols), dM/dt is negative: the magnetization 
decreases as population is moved from m = 10 to 9. After the 
radiation is turned off (open symbols), dM/dt becomes posi- 
tive as the population returns to the ground state at the ~5-/US 
phonon bottleneck time scale. 

When we raise the temperature to 3.3 K, the dynamics be- 
come richer, as shown in Fig. 3b. Immediately after the radia- 
tion is turned on (0.145 /j,s), the lineshape looks approximately 
symmetric, similar to the 2.1 K data. However, at later times 
a shoulder develops on the high-field side of the curve. After 
the radiation is turned off, dM/dt in the vicinity of the shoulder 
remains negative with a minimum at '^2650 Oe. Thus, the mag- 
netization continues to decrease after the radiation is turned off. 
The position of this minimum corresponds to the field at which 
resonant tunneling occurs between levels in opposite wells for 
Feg at the crystal's orientation. 

These results can be understood in terms of a thermally 
assisted resonant tunneling process driven by the radia- 
tion/phonon bottleneck, as follows (Fig. 4 insets). The radiation 
promotes some population to the m = 9 state, where it stays for 
some time because of the phonon bottleneck. If the tempera- 
ture is low (e.g. 2.1 K), there are insufficient thermal phonons 
to significantly populate any other level on the order of Tpb and 
the system decays back to equilibrium. At the higher tempera- 
ture (e.g. 3.3 K), thermal phonons are able to excite the system 
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Fig. 2: (Color online) Induced emf for 2-/is radiation pulses applied to 
the Feg sample at a) 2.1 K and b) 3.3 K at several values of magnetic 
field, as noted. We define t=0 to coincide with the time when the 
radiation is turned on. Data at many values of magnetic field were 
taken; only a small subset is shown for clarity of presentation. Each 
curve represents the average of 12800 oscilloscope traces in a) and 
16000 traces in b). The time between pulses was 20 ms. 



to higher levels (m = 8, 7, . . . ). Some population is then trans- 
ferred into the opposite well by tunneling. Even a small popu- 
lation change between wells can lead to a relatively large signal 
since Am is large ('-^10-20) for tunneling transitions. The role 
of the phonon bottleneck is essential to the dynamics: it slows 
down the transitions between the m = 9 and 10 states, which 
is primarily responsible for the 2.1-K results, while at higher 
temperatures (e.g. 3.3 K) thermally assisted tunneling effects 
can be observed on faster ('^100-ns) timescales because there 
is no bottleneck for transitions between higher-lying states. 

We can quantitatively model these dynamics through a 
master-equation approach [21-24] that includes radiative tran- 
sitions, phonon transitions and treats the population of phonons 
resonant with the lO-to-9 transition as a dynamical variable. 
We work in the spin's energy eigenbasis, which incorporates 
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Fig. 3: (Color online) Induced emf as a function of magnetic field for 
several points in time after the onset of a 2- /is radiation pulse, as noted, 
obtained from the data in Fig. 2 and similar data not shown, at a) 2. 1 
K and b) 3.3 K, averaged over a 100-ns time window. 

tunneling effects automatically through the fact that the eigen- 
states of Eq. 1 are superpositions of rn levels. We neglect off- 
diagonal elements in the density matrix, which is a good ap- 
proximation except extremely close to a tunneling resonance 
(when two levels are within a natural linewidth of each other). 
The rate of change of the population of the ith energy eigenstate 
is given by: 



(2) 

The first term represents transitions from the ith level to the oth- 
ers and the second term represents transitions to the ith level. 
The populations are normalized so that = 1. The pho- 

i 

ton (wij) and phonon (7,^) transition rates are calculated using 
golden-rule expressions [18,22]: 
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Fig. 4: (Color online) Results of simulations. The calculated values 
of dM/dt as a function of field at various times after the onset of a 
radiation pulse. Symbols are the same as in Fig. 3. The values of 
dM/dt have been scaled by experimentally determined parameters to 
have the units of induced emf ipW). The insets show schematically the 
dominant level transitions determined by the simulations. Red (green) 
arrows represent primarily photon- (phonon-) driven transitions. 
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2h^ 
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\{mS±.S.Mj)\\e,-e,fN^I', 



(3) 



where S± are the standard spin raising and lowering operators, 
Ei is the energy of level p is the mass density, and is 
the transverse speed of sound. F{lo) is the resonance lineshape 
function, which we take to be a Gaussian of width a. We do 
not consider collective spin-phonon interactions [25]. The rates 
Wij are appreciable only during the radiation pulse [26] . For the 



phonon transition rates in Eq. 3, N^^ 



(+1) for absorption 



(per spin) resonant with the i ^ j transition. For most of the 
phonon transitions we set Nfj^ to its thermal equilibrium value 
(•g(ei-ej)/fcBT _ 1)^1. For transitions between the m = 10 (jg)) 
andm = 9 (|e}) states, we treat the phonon number as a variable 
ripb, the phonon-bottleneck number: 



pb 



dt 



P. 



[npb + 1) - 



Tpb 



(emission) of phonons, where is the number of phonons 



Here the first term on the right-hand side is the rate of phonon 
absorption, the second term, the rate of phonon emission, in- 
cluding spontaneous emission, and the last term, the rate at 
which the phonons decay toward their equilibrium population. 
We numerically solve Eqs. 2 and 4 for the populations Pi and 
Upb as functions of time, with thermal equilibrium values used 
as initial conditions. The rate of magnetization change is calcu- 

lated using = g (,| 5? , ^ where 5 • ^ is the 

projection of the spin operator along the external-field (mea- 
surement) direction. Finally, the results are convoluted with 
a Gaussian of width 200 Oe to account for inhomogeneous 
broadening of the resonant tunneling features. The results of 
these simulations are shown in Fig. 4. 

In performing the simulations, we used the known Hamil- 
tonian and the spectroscopically determined orientation of the 
sample. The position, width and depth of the main (pho- 
ton/bottleneck) peak are respectively determined by parameters 
D, a and Hx\ these we set to 0.290 K, 650 Oe and 1.0 Oe, re- 
spectively, each in good agreement with its independently mea- 
sured value {D and a determined from ESR spectra.). We set 
the phonon bottleneck time to the measured value of Tph = 5 /zs 
(Fig. 1). We treat the temperature T as a constant on the 
timescale probed here. We set p = 1920 kg/m^^ [27]The time de- 
pendence of the magnetization in our simulations is extremely 
sensitive to the value of Cs (because of the factor in Eq. 3). 
We found = 670 m/s to give the best agreement between 
the simulations and our 3.3-K data. This value is 16% smaller 
than the value extracted from specific heat measurements [28]; 
the discrepancy could arise from the fact that specific heat 
measures the average speed of sound, not the transverse com- 
ponent. It has been proposed that spin-phonon transitions in 
which Am = 2 could provide an additional relaxation mecha- 
nism and increase the relaxation rate [23,24]. We were unable 
to adequately simulate our data using such a model. Using the 
above parameters, we calculated the system's relaxation rate in 
the absence of radiation near zero magnetic field at several tem- 
peratures and obtained results consistent with those found in ac 
susceptibility measurements on the same material [29]. From 
Eq. 3 and the fit value of c^, we deduce a zero-temperature (i.e. 
spontaneous emission) lifetime of the m = 9 excited state of 
Ti ^ 40 ns. The simulations also allow us to deduce the dom- 
inant level transitions, shown in the insets to Fig. 4, with the 
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dominant tunneling transition occurring between levels m = 6 
and -7. 

Our results are in contrast to some recent work by another 
group studying the same material. Petukhov et al. [10] have 
suggested that a phonon bottleneck is responsible for the >1- 
ms radiation-induced dynamics in Peg. The time scale of those 
observations, however, corresponds to the heating effect char- 
acterized in [6, 8]. Very recently, Bahr et al. [12] used a pump- 
probe technique to study the magnetization dynamics of Fes 
exposed to pulsed radiation. They inferred from their data that 
Ti ~ 2/is for the m = 9 state, nearly two orders of magnitude 
slower than the value determined in this study, but close to our 
value of Tpb. We found it impossible to simulate the tunneUng 
results (Fig. 3b) with such a long Ti, essentially because the 
thermal phonon transition rates must be slower than 1/Ti, pre- 
venting thermally assisted tunneUng dynamics from occurring 
on the timescale of ^100 ns when Ti ^ 2/is, even if Am = 2 
processes are included. The discrepancy between the two sets 
of experiments is the subject of ongoing investigation. 

Our simulations capture all of the qualitative features of our 
data and give very good quantitative agreement with the 3.3- 
K results and approximate agreement with the 2.1-K data. We 
note that the peak in the 2.1-K data is shifted slightly toward 
lower field relative to its position at 3.3 K, an effect that has 
been ascribed to dipole fields [30]. While the simulations rea- 
sonably reproduce the 2.1-K data at fields below '^1000 Oe, 
they show a larger signal and stronger time dependence near 
the minimum than the actual data. We attribute this to the fact 
that when the sample is on resonance with the cavity (i.e. near 
the peak), its absorption effectively reduces the cavity Q, low- 
ering the value of Hi and therefore the photon absorption rate. 
This effect is more pronounced at lower temperature when the 
ground-state population is larger. It will also affect the time 
dependence as the populations change with photon absorption. 

The line width a of the photon-absorption resonances rep- 
resents the effect of inhomogeneous broadening. Our observa- 
tion of tunneUng constrains the physical nature of this broaden- 
ing. The 3.3-K results require a coincidence between two reso- 
nances: the photon resonance where the photon energy matches 
that of the lO-to-9 transition and the tunneUng resonance in 
which levels in opposite wells align. This obviates the pos- 
sibility that the broadening a could arise primarily from some 
quasi-static inhomogeneous dipole or hyperfine field since such 
a random field would shift both resonances equally and could 
not result in their overlap. On the other hand, our results are 
consistent with a model in which there is a distribution in val- 
ues of the anisotropy parameter D [30], wherein only a fraction 
of the molecules have a value of D that satisfies the coincidence 
conditions. 

In summary, we have measured the fast magnetization dy- 
namics of the single-molecule magnet Fes subject to short 
pulses of intense microwave radiation. We observe the develop- 
ment of a phonon bottleneck with a decay time of ^5 /is. The 
bottleneck can drive resonant tunneling between excited states 
when thermal phonons are available to populate these states, 
i.e. at sufficiently high temperature. At low temperatures, the 
dynamics are restricted to the the two states resonant with the 



radiation and dominated by the phonon-bottleneck time scale. 

Our detailed numerical simulations are in good agreement with 
the data and imply an excited state lifetime of ~40 ns. 
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